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1.  Introduction  > 

'  Irrthis  paperr  we-discuss  recent  results  on  the  nonlinear  stability  of  discontinuous 
steady  states  of  a  mod^f  initial-boundary  value  problem  in  one  space  dimension  for 
incompressible,  isothermal  shear  flow  of  a  non-Newtonian  fluid  between  parallel  plates 
located  at  x  =  ±1,  and  driven  by  a  constant  pressure  gradient.  The  non-Newtonian 
contribution  to  the  shear  stress  is  assumed  to  satisfy  a  simple  differential  constitutive 
law.  The  key  feature  is  a  non-monotone  relation  between  the  total  steady  shear  stress 
and  steady  shear  strain  rate  that  results  in  steady  states  having,  in  general,  discontinu¬ 
ities  in  the  strain  rate.  We  explain  why  every  solution  tends  to  a  steady  state  as  t  — *■  bo, 
and  we  identify  steady  states  that  are  stable;  more  details  and  proofs  will  be  presented 
in  [8].  .  ]  f 

We  study  the  system  p  ( - 


(1.1) 


Vt  =  S : 


r  > 


S  :=  T  +  /x, 


T  :=  cr  +  v 


X  > 


(1.2)  <7t  +  CT  =  g{vs)  , 

on  [0, 1]  x  [0,  oo),  with  /  a  fixed  positive  constant.  We  impose  the  boundary  conditions 

(1.3)  5(0, f)  =  0  and  u(l,t)  =  0,  t  >  0, 
and  the  initial  conditions 

(1.4)  u(x,0)  =  v0(z),  cr(z,  0)  =  <70(x),  0<z<l; 

accordingly,  5(x,  0)  =  5o(x)  :=  <7q(x)  +  uox(x)+/x.  The  function  g  :  IR  — »  1R  is  assumed 
to  be  smooth,  odd,  and  £g(£)  >  0,  ^  ^  0.  In  the  context  of  shear  flow,  v,  the  velocity 
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of  the  fluid  in  the  channel,  and  T,  the  shear  stress,  axe  connected  through  the  balance 
of  linear  momentum  (1.1).  The  shear  stress  T  is  decomposed  into  a  non-Newtonian 
contribution  a,  evolving  in  accordance  with  the  simple  differential  constitutive  law  (1.2), 
and  a  viscous  contribution  vx.  The  coefficients  of  density  and  Newtonian  viscosity  are 
taken  as  1,  without  loss  of  generality.  The  flow  is  assumed  to  be  symmetric  about  the 
centerline  of  the  channel.  Symmetry  dictates  the  following  compatibility  restrictions  on 
the  initial  data: 

(1.5)  v0(l)  =  0,  So(0)  =  0,  and  a0(  0)  =  0; 

they  imply  that  <x(0,i)  =  ux(0,t)  =  0,  and  symmetry  is  preserved  for  all  time. 

The  system  (1.1)-(1.4)  admits  steady  state  solutions  (u(x),ct(x))  satisfying 


(1.6)  S  :=  g(yx)  +  vx  +  fx  =  0,  a  =  g(vx) 

on  the  interval  [0, 1].  In  case  the  function  «;(£)  :=  g(£)  +  f  is  not  monotone,  there 
may  be  multiple  values  of  vx(x)  that  satisfy  (1.6)  for  some  x's,  thus  leading  to  steady 
velocity  profiles  with  jumps  in  the  steady  velocity  gradient  vx.  Our  objective  is  to 
study  the  stability  of  such  steady  velocity  profiles;  we  also  study  well-posedness  and  the 
convergence  of  solutions  of  (i.l)-(1.4)  to  steady  states  a s  t  oo. 


w£) 


a  representative  w(£)  is  shown  in  Fig.  1;  in  the  figure  m  and  M  stand  for  the  levels  of 
the  bottom  and  top  of  the  loop,  respectively.  Our  results  and  techniques  can  be  easily 
generalized  to  cover  the  case  when  u>(£)  has  a  finite  number  of  loops.  Steady  state 
velocity  profiles  are  constructed  as  follows:  First  solve  w(u(x))  =  fx  for  each  x  6  [0, 1], 
where  u  =  — vx.  This  equation  admits  a  unique  solution  for  0  <  fx  <  m  or  fx  >  M , 
and  three  solutions  for  m  <  fx  <  M;  let  u(x),  0  <  x  <  1,  be  a  solution.  Setting 

(1.7)  v(x)  =  J  u(y)dy ,  V(x)  =  p(-u(x)) , 

then  (u(x),  &(x))  satisfy  (1.6)  and  (1.3)  for  a.e.  x  €  [0, 1]  and  give  rise  to  a  steady  state. 
Clearly,  if  /  <  m  there  is  a  unique  smooth  steady  state;  if  m  <  /  <  M,  there  is  a 


"or 


sa 


□ 

□ 


-n/ _ 

ty  Codes 

•ind/or 

lal 


2 


unique  smooth  velocity  profile  and  a  multitude  of  profiles  with  kinks;  finally,  if  /  >  M, 
all  steady  state  velocity  profiles  have  kinks.  An  example  of  a  velocity  profile  with  kinks 
is  shown  in  Fig.  2. 

Problem  (1.1)-(1.4)  captures  certain  key  features  of  a  class  of  viscoelastic  mod¬ 
els  that  have  been  proposed  to  explain  the  occurence  of  “spurt”  phenomena  in  non- 
Newtonian  flows.  Specifically,  for  a  particular  choice  of  the  function  g  in  (1.2),  the 
system  under  study  has  the  same  steady  states  as  the  more  realistic  systems  studied  in 
[6]  and  [7];  the  latter,  derived  from  a  three-dimensional  setting  that  is  restricted  to  one¬ 
dimensional  shearing  motions,  produce  non-monotone  steady  shear  stress  vs.  strain-rate 
relations  of  the  type  shown  in  Fig.  2.  The  phenomenon  of  spurt  was  apparently  first 
observed  by  Vinogradov  et  al.  [13]  in  the  flow  of  highly  elastic  and  very  viscous  non- 
Newtonian  fluids  through  capillaries  or  slit-dies.  It  is  associated  with  a  sudden  increase 
in  the  volumetric  flow  rate  occuring  at  a  critical  stress  that  appears  to  be  independent 
of  the  molecular  weight.  It  has  been  proposed  by  Hunter  and  Slemrod  [5],  using  tech¬ 
niques  of  conservation  laws,  and  more  recently  by  Malkus,  Nohel,  and  Plohr  [6]  and  [7], 
using  numerical  simulation  and  extensive  analysis  of  suitable  approximating  dynamic 
problems  (motivating  the  present  work),  that  spurt  phenomena  may  be  explained  by 
differential  constitutive  laws  that  lead  to  a  non-monotone  relation  of  the  total  steady 
shear  stress  versus  the  steady  shear  strain-rate.  In  this  framework,  the  increase  of  the 
volumetric  flow  rate  corresponds  to  jumps  in  the  strain  rate  when  the  driving  pressure 
gradient  exceeds  a  critical  value.  We  conjecture  that  our  stability  result  discussed  in 
Sec.  3  below  can  be  extended  to  these  more  complex  problems. 

2.  Preliminaries 

In  this  section,  we  discuss  preliminary  results  that  are  essential  for  presenting  the 
stability  result;  further  details  and  proofs  can  be  found  in  [8]. 

A.  Well-Posedness. 

We  use  abstract  techniques  of  Henry  [4]  to  study  global  existence  of  classical  so¬ 
lutions  for  smooth  initial  data  of  arbitrary  size ,  and  also  existence  of  almost  classical, 
strong  solutions  with  discontinuities  in  the  initial  velocity  gradient  and  in  the  stress 
components.  The  latter  result  allows  one  to  prescribe  discontinuous  initial  data  of  the 
same  type  as  the  discontinuous  steady  states  studied  in  this  paper.  Ex;  t..ice  results  of 
this  type  are  established  in  [8]  for  a  general  class  of  problems  that  serw  5.3  models  for 
shearing  flows  of  non-Newtonian  fluids;  the  total  stress  is  decomposed  into  a  Newtonian 
contribution  and  a  finite  number  of  stress  relaxation  components,  viewed  as  internal 
variables  that  evolve  in  accordance  with  differential  constitutive  laws  frequently  used 
by  rheologists  (for  discussion,  formulation  and  results,  see  [11],  [7],  also  the  Appendix 
in  [8]).  Existence  of  classical  solutions  may  also  be  obtained  by  using  an  approach  based 
on  the  Leray  -  Schauder  fixed  point  theorem  (cf.  Tzavaras  [12]  for  existence  results  for 
a  related  system).  Other  existence  results  were  obtained  by  Guillope  and  Saut  [2],  and 
for  models  in  more  than  one  space  dimension  in  [3]. 

As  a  consequence  of  the  general  theory,  one  obtains  two  global  existence  results 
(see  Theorems  3.1,  3.2,  3.5,  and  Corollary  3.4  in  [8]): 

(i.)  the  existence  of  a  unique  classical  solution  (v(x,t),o-(x,t))  of  (1.1)-  (1.5)  on  [0, 1]  x 

[0,oo)  for  initial  data  (uq(^),  <7o(z))i  n°t  restricted  in  size,  that  satisfy:  Sq(x )  := 
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Vox(x)+crQ(x)+fx  €  Hs[ 0, 1]  for  some  s  >  3/2,  with  5o(0)  =  0,  u0(l)  =  50r(l)  =  0> 
and  (7 o  G  C1  [0, 1],  where  H 3  denotes  the  usual  interpolation  space. 

(ii.)  the  existence  and  uniqueness  of  a  strong,  “semi-classical”  solution  of  (1.1)  -(1.5), 
obtained  by  a  different  choice  of  function  spaces,  for  initial  data  (vQ(x),  <r0(x))  that 
satisfy:  5o(x)  €  FT^O,  1]  with  5o(0)  =  0,  i/0(l)  =  50r(l)  =  0,  and  aQ  €  X°°[0, 1]. 
Result  (ii.)  yields  solutions  in  which  a  and  vx  may  be  discontinuous  in  x,  but  Sz  and  vt 
are  continuous,  and  a  is  C1  as  a  function  of  t  for  every  x.  Thus  all  derivatives  appearing 
in  the  system  may  be  interpreted  in  a  classical  sense  as  long  as  the  equation  is  kept  in 
conservation  form.  A  result  of  this  type  was  obtained  by  Pego  in  [10]  for  a  different 
problem  by  a  similar  argument. 

B.  A  Priori  Bounds  and  Invariant  Sets. 

To  discuss  global  boundedness  of  solutions,  let  <r,t)  be  a  classical  solution  on  an 
arbitrary  time-interval,  and  note  that  system  (1.1)  -(1.5)  is  endowed  with  the  differential 
energy  identity 


(2.1) 


—  {1/2  J  v*dx  +  J  [W(vx)  +  x  f  vx]dx} 

+  [  \vt  +  vlt\dx  =  0  • 

Jo 


The  function  W(£)  :=  w(Qd(  plays  the  role  of  a  stored  energy  function;  by  the 
assumption  on  g,  W  is  not  convex.  This  fact  is  the  main  obstacle  in  the  analysis  of 
stability. 

{i.) Boundedness  of  S.  Since  £g(£)  >  0,  it  follows  that  g(Qd(  >  0  for  £g!R,  and 
W( £)  satisfies  the  lower  bound 

(2.2)  W(£)  +  fx$>l/4e-f2,  eelR,  0  <  x  <  1 . 


Standard  energy  estimates  based  on  (2.1)  and  (2.2)  coupled  with  integration  of  (1.1) 
with  respect  to  x  yield  a  global  a  priori  bound  for  S: 


(2.3)  |S(x,f)|  <  C  0  <  x  <  1,  0  <  t  <  oo , 

where  C  is  a  constant  depending  only  on  data  but  not  t. 

(ii.) Invariant  Sets  for  a  Related  ODE.  Control  of  S  enables  us  to  take  advantage 
of  the  special  structure  of  Eq.  (1.2)  and  determine  suitable  invariant  regions.  For  this 
purpose,  it  is  convenient  to  introduce  the  quantity  s  :=  a  +  f  x.  Then,  Eqs.  (1.2),  (1.3) 
readily  imply  that  s  satisfies 


(2.4)  st  +  s  +  g(s-S)  =  fx. 

For  a  fixed  x,  it  is  convenient  to  view  Eq.  (2.4)  as  an  ODE  with  forcing  term  5(x,«). 
Also,  observe  that  at  a  steady  state  (a,vx),  one  has  5  =  0,  and  consequently,  s  =  —vz 
is  an  equilibrium  solution  of  (2.4)  (with  5  =  0).  If  5  s  0  in  (2.4),  the  hypothesis 
ronceming  g  implies  that  the  ODE  admits  positively  invariant  intervals  for  each  fixed 
x.  We  sketch  how  this  property  is  preserved  in  the  presence  of  a  priori  control  oi  S  as 
provided  by  (2.3);  more  delicate  bounds  axe  essential  in  the  proof  of  stability  in  Sec  .3. 
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To  fix  ideas,  let  to  >  0  be  given,  and  assume  that 

(2.5)  |S(s,i)|  </>,  0<ar<l,0<f<£o, 

for  some  p  >  0.  For  x  fixed  in  [0, 1],  we  use  the  notation  S(t)  :=  5(x, f)  and  conveniently 
rewrite  (2.4)  as 

(2.6)  st  +w(s  -  S(t))  =  fx  —  S(t) . 

We  state  the  following  result  on  invariant  intervals;  its  proof  is  obvious. 

Proposition  2.1.  Let  S  satisfy  the  uniform  bound  (2.5)  for  0  <  t  <  to.  For  x  fixed , 
0  <  x  <  1,  assume  there  exist  s_,  s+  such  that  s_  <  s+  and 


(2.7) 

w($~  —  A)  <  /  x  —  A  , 

|A|<P 

(2.8) 

w(s+  —  A)  >  /  x  —  A  , 

Then  the  compact  interval  [s_,  s+]  is  positively  invariant  for  the  ODE  (2.6)  on  the  time 
interval  0  <  t  <  to- 

Invariant  intervals  axe  generated  by  solution  sets  of  the  inequalities  (2.7)  and  (2.8) 
as  functions  of  p  and  x.  In  particular,  since  lim  u>(£)  =  ±oo,  given  any  x  and  p ,  one 

*±oo 

easily  determines  s0+  large,  positive  and  so-  large,  negative  such  that  if  s_  <  so-  and 
s+  >  so+,  then  s_  and  s+  satisfy  (2.7)  and  (2.8),  respectively,  and  the  compact  interval 
[s_,s+]  is  positively  invariant  for  the  ODE  (2.6). 


W  (5) 


Fig.  3:  Invariant  Intervals. 
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More  discriminating  choices  of  invariant  intervals  occur  if  one  restricts  attention  to 
small  values  of  p\  the  analysis  becomes  more  delicate.  For  a  function  u>(£)  with  a  single 
loop,  the  most  interesting  case  arises  when  /  x  —  p,  f  x  and  /  x  +  p  each  intersects  the 
graph  of  w(£)  at  three  distinct  points.  Referring  to  Fig.  3,  the  abscissae  of  the  points 
of  intersection  axe  denoted  by  (a_, /?_,7_),  (a o, $>,70)  (ar+,/?+,7+),  respectively. 

It  turns  out  that  for  x  fixed  and  p  small  enough,  there  are  discriminating  invariant 
intervals  of  the  type  shown  in  Fig.  3.  However,  in  contrast  to  the  large  invariant  intervals 
discussed  in  the  previous  paragraph,  the  more  discriminating  ones  degenerate  as  we 
approach  the  top  or  bottom  of  the  loop  (when  x  varies). 

For  the  stability  of  discontinuous  steady  states  in  Sec.  3.,  it  is  crucial  to  construct 
compact  invariant  intervals  that  are  of  uniform  length  (see  Corollary  2.2  in  [8]).  The 
latter  is  accomplished  by  taking  p  sufficiently  small  and  by  avoiding  the  top  and  bottom 
of  the  loop  in  Fig.  3.  Of  specific  interest  is  the  situation  in  which  s(x)  is  a  piecewise 
smooth  solution  of 

(2.9)  w(s(x))  =  fx, 

defined  on  [0, 1]  and  admitting  jump  discontinuities  at  a  finite  number  of  points  x1}...,xn 
in  [0, 1].  Recall  that  s(x)  is  a  steady  solution  of  the  ODE  (2.6)  corresponding  to  the 
steady  state  (<t,vx).  In  addition,  suppose  that  s(x)  takes  values  in  the  monotone  in¬ 
creasing  parts  of  the  curve  io(£)  and  that  it  avoids  jumping  at  the  top  or  bottom  of  the 
loop,  i.e., 

(2.10)  w'(s(x))  >  c0  >  0 ,  x  6  [0, 1]  \  {xi,...,xn}  , 

for  some  constant  c0.  A  delicate  construction  in  [8]  yields  compact,  positively  invariant 
intervals  of  (2.6)  of  uniform  length,  centered  around  s(x)  at  each  x  €  [0, 1] \  {xj, ..., x„}. 

(ni.)  Boundedness  of  a  and  vz.  As  an  easy  application  of  Sec.  2  (ii),  choose  a 
compact  interval  [s_,  s+]  that  is  positively  invariant  for  (2.6)  and  valid  for  all  x  €  [0, 1]. 
By  virtue  of  the  global  bound  (2.3)  satisfied  by  S(x,<),  we  conclude  that 

(2.11)  \s(x,t)\<C,  0<x<l,t>0 

which,  in  turn,  using  (1.3)  and  (2.11),  implies 

(2.12)  Mx,t)|<C,  0<x<l,  <>0, 

for  some  constant  C  depending  only  on  the  data.  The  definition  of  s  also  implies  that 
<7  is  uniformly  bounded. 

(iv .)  Convergence  to  steady  states.  Let  (v(x,t),  a(x,t))  be  a  classical  solution  of 
(1.1)  -(1.5)  defined  on  [0, 1]  x  [0,  00].  We  discuss  the  behavior  of  this  solution  as  t  — ►  00. 

The  first  result  indicates  that  S  =  <7  +  vI  +  /x  converges  to  its  equilibrium  value. 
Proposition  2.2.  Under  the  assumptions  of  the  existence  results , 

(2.13)  lim  5(x,t)  =  0 , 

'  '  *  — *00 

uniformly  for  x  €  (0, 1]. 
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The  proof  is  a  consequence  of  Sobolev  embedding  applied  to  the  following  a  priori 
estimates  that  axe  derived  from  the  system  (1.1)— (1.4)  by  standard  techniques: 


(2.14) 

[°°  f  Sfdxdr  <  C  , 
Jo  Jo 

(2.15) 

f  [  S2dxdr  <  C , 
Jo  Jo 

(2.16) 

f  Sl(x,t)dx  <  C,  0  <  t  < 
Jo 

where  C  is  a  positive  constant  depending  only  on  data. 

Use  of  (2.13)  enables  us  to  identify  the  limiting  behavior  of  solutions  of  (2.4)  as 
t  — >  oo.  The  following  result  is  analogous  to  Lemma  5.5  in  Pego  [10];  its  elementary 
proof  is  given  in  Lemma  4.2  of  [8]. 

Proposition  2.3.  Let  s(x,«)  €  C1[0,oo)  be  the  solution  of  (2.10),  where  S(x,»)  is 
continuous  and  satisfies  (2.13),  0  <  x  <  1.  Then  s(x,  •)  converges  to  s°°(x)  as  t  — *  oo 
and  s°°(x)  satisfies 

(2.17)  s°°(x)  -f-  y(s°°(x))  —  f  x  ,  0  <  x  <  1 . 

In  view  of  the  shape  of  te(<f)  =  £  +  g(£),  equation  (2.17)  has  one  solution  for  0  <  fx  <  m 
or  f  x  >  M  and  three  solutions  for  m  <  f  x  <  M . 

Let  (u(x,  t),  cr(x,  t))  be  a  classical  solution  of  (1.1)-(1.4)  on  [0, 1]  x  [0,  oo).  Recalling 
the  definition  of  s,  Proposition  2.3  implies 

(2.18)  cr°°(x)  =  lii^  cr(x,  t )  =  s°°(x)  -  /  x  . 

Also,  combining  (1.1),  (2.13)  and  (2.18)  yields 

(2.19)  v£°(x)  :=  lim  vx(x,t)  =  lim  (S(x,  t )  -  s(x,  <))  =  - s°°(x ) , 

N  t— ♦  oo  t— *oo 

and 

(2.20)  5°°(x)  =  t£°(x)  +  cr°°(x)  +  fx  =  0. 

Finally,  noting  that 

(2.21)  0  =  J  vx(x,  t)dx  , 

u°°(x)  is  Lipschitz  continuous  and  satisfies 

(2.22)  u°°(x)  :=  lim  v(x,  t)  =  [ 

‘—o°  Jx 
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We  conclude  that  any  solution  of  (1.1)-(1.4)  converges  to  one  of  the  steady  states. 
If  0  <  /  <  m,  then  there  is  a  unique  smooth  steady  state  which  is  the  asymptotic  limit 
of  any  solution.  However,  if  m  <  /,  then  there  are  multiple  steady  states  and  thus  a 
multitude  of  possible  asymptotic  limits.  In  Sec.  3,  we  identify  stable  steady  states.  Also 
note  from  (2.20)  that  in  a  discontinuous  steady  state,  the  discontinuities  in  a  and  vx 
cancel. 

Observe  that  in  case  u>(£)  is  monotone  the  above  arguments  yield  that  every  solution 
converges  to  the  unique  steady  state.  Moreover,  the  above  results  can  be  routinely 
generalized  to  the  case  that  the  function  u>(£)  has  multiple  loops  but  the  graph  of  w 
has  no  horizontal  segments. 

3.  Stability  of  Steady  States 

The  purpose  is  to  study  the  stability  of  velocity  profiles  with  kinks.  To  fix  ideas, 
let  (u(x),<t(x))  be  a  steady  state  of  (1.1) — (1.3)  such  that  u(x)  has  a  finite  number  of 
kinks  located  at  the  points  x\, . . . ,  xn  in  (0, 1);  accordingly,  vx(x)  and  er(x)  have  a  finite 
number  of  jump  discontinuities  at  the  same  points.  Recall  that,  if  we  set  ll(x)  =  —vx(x), 

(3.1)  w(u(x))  =  fx,  x  €  [0, 1],  x  ^xi,...,xn 
and  <r(x)  =  g(—u(x)). 

Given  smooth  initial  data  (u0(x),  <To(x)),  there  is  a  unique  smooth  solution  (u(x,t), 
a(x,t))  of  (1.1)-(1.4).  As  t  — »  oo,  the  solution  converges  to  one  of  the  steady  states, 
not  a-priori  identifiable.  We  now  restrict  attention  to  initial  data  that  are  close  to 
(u(x),  a(x)),  except  on  the  union  U  of  small  subintervals  centered  around  the  points 
Xi,...,x„.  U  can  be  thought  of  as  the  location  of  transition  layers  separating  the 
smooth  brandies  of  the  steady  state.  Roughly  speaking,  it  turns  out  that  the  steady 
state  is  “asymptotically  stable”  under  smooth  perturbations  that  are  close  in  energy, 
provided  (u(x),  <f(x))  takes  values  in  the  monotone  increasing  parts  of  ir(£);  the  stable 
solutions  are  local  minimizers  of  an  associated  energy  functional  (see  (3.8)  below).  The 
interesting  problem  of  finding  the  domain  of  attraction  of  a  stable  steady  solution  ap¬ 
pears  to  be  a  difficult  task.  Our  main  result  is: 

Theorem  3.1.  Let  (u(x),a(x))  be  a  steady  state  solution  as  described  above  and  satis¬ 
fying 

(3.2)  ty'(us(x))  >  Co  >  0,  x  G  [0,1],  x^xx,...,xn 

for  some  positive  constant  Co ■  If  the  measure  ofU  is  sufficiently  small,  there  is  a  positive 
constant  Sq  depending  onU  such  that,  if  8  <  <50,  then  for  any  initial  data  (vo(x),0'q(x)) 
satisfying 

(3.3)  sup  j50(x)|  <  8, 

0<r  <1 

r  1  2 

(3.4)  J  vf(x,0 )dx  < -82 


8 


and 

(3.5)  |uor(z)  ~vx(x)\  <  8  ,  i6[0,l]\W 

the  corresponding  solution  (u(x,f),a(x,<))  approaches  the  steady  state  (u(x),  ct(x))  as 
t  — ►  oo,  in  the  sense, 


(3.6) 


V*(x,f)  ->Vx(x), 


(3.7)  cr(x,t)  -»■  <r(x), 

/or  a//  x  6  [0, 1]  \  U . 

The  above  result  is  similax,  in  spirit  and  approach,  to  the  analysis  of  Andrews  and 
Ball  [1],  and  particularly  to  stability  results  established  by  Pego  [10]  for  motions  of 
one-dimensional  viscoelastic  materials  of  rate  type  with  a  non-monotonic  stress-strain 
relation.  Current  work  of  Novick-Cohen  and  Pego  [9]  on  spinodal  decomposition  involve 
a  similax  stability  analysis. 

Because  the  argument  is  lengthy  and  delicate,  we  can  only  indicate  the  main  idea 
of  the  proof  for  the  case  that  u(x)  =  — ur(x)  has  one  single  jump  discontinuity  located 
at  xo,  m  <  fx0  <  M,  and  for  U  =  (x<j  —  £,x0  +  e)  for  some  small  e.  Minor  modifica¬ 
tions  axe  needed  to  account  for  the  general  case.  Technical  details  can  be  found  in  [8], 
Theorem  5.1.  The  proof  is  based  on  exploiting  the  energy  identity  (2.1),  which,  upon 
setting  u(x,t)  —  —vx(x,t)  and  integrating  with  respect  to  t ,  yields  the  inequality 


(3.8) 


^  Jo  +  JQ  [(W(u(x,t))  -  xfu(x,t))  -  $(x)]dx 

[  Uj(x,0)dx  +  f  [(W(u0(x))  —  xfuo(x)  —  $(x)]dx; 

^  Jo  Jo 


note  that  the  integral  of  the  function  $  has  been  subtracted  from  both  sides  of  (3.8). 
The  function  $(x)  is  associated  with  the  particular  choice  of  the  discontinuous  steady 
state  u(x),  the  stability  of  which  is  being  tested.  The  function  $  identifies  a  basin  of 
attraction  of  the  state  u(x).  Roughly  speaking,  the  goal  is  to  find  $  so  that  the  second 
integral  on  the  left  side  of  (3.8)  is  positive,  and  at  the  same  time  the  right  side  can 
be  made  sufficiently  small.  The  construction  of  <5  is  delicate  because  W  in  (3.8)  is  not 
convex,  and  because  the  double- well  potential  W(u)  —  xfu  depends  explicitly  on  x. 
Note  that  for  each  fixed  x,  the  function  u(x)  is  the  horizontal  coordinate  of  the  bottom 
of  one  of  the  wells;  the  left  well  if  x  <  xo,  the  right  well  if  x  >  xo-  To  insure  that  the 
construction  of  $  produces  the  property  desired,  this  paxt  of  the  analysis  makes  crucial 
use  of  invariant  intervals  of  the  ODE  (2.6)  that  are  of  uniform  length  as  discussed  in 
Sec.  2(ii)  above. 
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